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ABSTRACT 

We propose a new model for the dark matter halo of the Milky Way that fits the properties of the stellar stream 
associated with the Sagittarius dwarf galaxy. Our dark halo is oblate with q^ — 0.9 for r < 10 kpc, and follows 
the Law & Majewski (2010) model at large radii, the transition between these regimes occurring at ^ 30 kpc. The 
outer halo can be made mildly triaxial, with minor-to-major axis ratio (c/a)^ = 0.8 and intermediate-to-major 
axis ratio [h/a)^ — 0.9, if the effect of the Large Magellanic Cloud is taken into account. Therefore this new 
model takes into account the flattening induced by the presence of the Galactic disk and is also more consistent 
with cosmological expectations. 

Subject headings: galaxies: dwarf- galaxies: interactions - Local Group - Galaxy: halo 



1. INTRODUCTION 

The stellar stream associated with the Sagittarius (Sgr) dwarf 
galaxy has been extensively used to probe the mass distribution 
of the Milky Way (MW), particularly its dark halo. Despite 
many attempts, there is currently no fully satisfactory model of 
its shape based the dynamics of the stream. Extrem e oblate con- 
figurations have been ruled out (llbataet al.ll200ll) . while the tilt 
of the orbit al plane has been shown to req uire a mildly oblate 
halo by Joh nston. Law & Maiewskil ( l2005h . On t he other hand , 
the line of sight velo cities call fo r a prolate halo (Helmi 2004). 
This conundrum led dLaw & Majewski 2010, hereafter LMIO) 
to propose a triaxial dark halo for the MW, with axis ratios 
(claU > = 0.72 and (6/a)* = 0.99 (see also iDeg & Widrowl 
I2OI2I) . This model fits well all positional and kinematic infor- 
mation available. 

Although the halos assembled in ACDM are as a rule triaxial 



dJing & Sutol2002l;lAllgood et al.l2006tlSchneider. Frenk & Cold 
1201 2*). the configuration proposed by LMIO is rare: the halo 
is close to oblate, with a much smaller c/a than predicted 
in cosmological simulations for MW-mass halos, which have 
lc/a)<i, = 0.9 ± 0.1 over th e relevant distance range 
dHavashi. Navarro & Springell2007b . and difficult to understand 
from a physical point of view (its minor axis points almost to- 
wards the Sun, while the intermediate axis is perpendicular to 
the Galactic disk). Furthermore, the presence of the disk is 
expected to lead to a change in the inner halo shape towards 
a more oblate configuration dBrvan et al.ll2013l) . Finally, the 
disk's stability is not naturally ensured in the LMIO poten- 
tial, as there are no tu be orbits around the intermediate axis 
dDebattistaetal.ll2013h . 

In this Letter we take a fresh look at determining the shape 
of the MW halo from the Sgr streams dynamics. We consider 
the possibility that the shape of the halo varies with distance 
from the Galactic cent er, as expected in the context of ACDM 
dVera-Ciro et al.ll201 lb . Evidence sugge sting a halo with non- 
constant axis ratios has been reported by Baneriee & Jog|(l201 ll) 
using the flaring of the Hi layer of the MW disk. We present a 
new model that takes into account the effect of a baryonic disk 
in Section 2. Because of the cosmological rareness of the LMIO 
model, in Sec. 3 we explore the possibility that the dynamics of 
the Sgr stream may be explained through the combined effect 
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of the Large Magellanic Could (LMC) and a less axisymmetric, 
but more triaxial, outer halo. In that section we show that these 
models provide equally good fits to the dynamics of the young 
Sgr streams as the LMIO potential, and that older wraps may 
be used to distinguish amongst them. We finalize with a brief 
summary in Sec. 4. 

2. INNER HALO: ACCOUNTING FOR THE EFFECT OF THE 
GALACTIC DISK ON THE HALO SHAPE 

Next we present the characteristics of our Galactic potential, 
which includes a halo whose shape by construction, is oblate in 
the center, and triaxial at large radii. We then show the results 
of orbital integrations in this potential aimed at reproducing the 
properties of the Sgr stream. 

2.1. Description of the potential 

We model the Galactic potential with three components: a 
disk, a spherical bulge and a dark matter halo. The disk 
and bulge follow, respectively a Miyamot o-Nagai distribution 
(Mriis k =1O"M0, a = 6.5 kpc, h = 0.26 kpc. lMivamoto & Nagail 
119751) . and a Hemquist spheroid (A/buigo = 3.4 x 10^° Mq, 
c = 0.7kpc,|Hernquist 1990). 

Based on the arguments presented in the Introduction, we 
seek a halo potential that satisfies: 

1 . It is axisymmetric in the inner parts. This will guarantee 
the stability of the disk, as well as account for the effects 
of the baryonic disk on the dark halo. 

2. It is triaxial in the outskirts, and follows the LMIO 
model. 

3. It has a smooth transition between these two regimes. 

We choose to model s uch a profile using a modi fication of the 
algorithm presented bv lVogelsberger et al.l(l2008l) . Consider the 
spherical potential 

*.M-«haioln(r'+d'). (1) 

The geometrical properties of the potential are encapsulated in 
thevariabler — (a;^ +j/^ + z^)^/^. A replacement that satisfies 
the above requirements is r — > r, with 
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Note that for r <^ Va r w ta, and similarly r ^ Va r ^ tt- 
If TA is therefore defined as a quantity with axial symmetry, a 
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potential of the form $s(r) will behave axisymmetrically for 
r ^Ta- The properties of the mass distribution for r :§> Tq are 
encoded in the quantity tt which we define the same way as 
LMIO. In summary. 
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The Sun is 



assumed to be located along the x axis, at a; = —Rq, and the z 
axis to point perpendicular to the disk. The resulting potential 

$haio(a;, j/,2;) = $s(r(x,y, z)), (6) 

is axisymmetric for r <^ra, and triaxial for r ':^ Va- 

Fig. [T] shows different slices of the resulting potential. Here 
we have chosen t he flattening of the axisymmetric part to be 
qz = 0.9 (as in e.g. lJohnston. Law & Maiewskil2005h . The axis 
ratios for the triaxial component (gi, (72, 93), as well as its tilt 
(/), are taken from the LMIO model. The circular velocity Whaio 
is set to ensure that VdrciRe — 8 kpc) = 225.2 km s^^. The 
transition radius Va = 30 kpc, is selected such that the region of 
dominance of the disk resides inside the axisymmetric part of 
the halo potential. However, the effective transition between the 
axisymmetric and triaxial regions, occurs at a smaller radius, 
X 10 kpc. 
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Fig. 1. — Dai'k halo potential isocontours on the plane z = (top left), 
y = (bottom left) and x = (top right). For reference we have included the 
positions and directions of motion for the Sun (circle), Sgr (squai'e) and the 
LMC (diamond). The bottom right panel shows the circular' velocity profile 
■"circ for the disk (dotted blue), bulge (dashed green) and halo (dash dotted 
red). The halo makes a transition from oblate to triaxial at r^ = 30 kpc. 



2.2. Generating the stream 

In what follows we work on the assumption that the orbit 
of the center of mass tr aces the arms of the st ream. Although 
this is not strictly true ('Eyre & Binnev"2009'), it represents a 
reasonable first approximation (Law & Majewski 2010). With 
this caveat, we proceed to integrate test particles in the com- 
posite potential described above. For each particle we generate 
a set of initial conditions consistent with the present day six 
dimensional phase-space coordinates of the Sgr dwarf galaxy. 
More specifically, we sample each observable from a Gaussian 
distribution, with its mean and variance taken from the litera- 
ture. T he position is assumed to be at (/, b) — (5.6°, —14.2°) 
JMaie wski et al. 200 3|), the h eliocentric distance d — 25 ± 2 
kpc ( Kunder & Chab over 2009), the line of sight velocity Vr = 
140 ± 2 km s~^ d lbataetal. 199 7) and the p roper motions 
inico sh^Hb) = (-2.4±0.2. 2.1±0.2) mas yr~^ (iDinescu et alJ 
I2OO5I) . Orbits are integrated forward and backward in time for 
2 Gyr, to generate the set of observables associated with the 
leading and trailing arm, respectively. 

For each integrated orbit, we take 10 samples of the form 
{x{ti) , v{ti)}YLi, where the times ti are randomly selected be- 
tween t ~{) and the maximum time of integration imax- imax is 
the time that it takes the orbit to complete one wrap in the sky, 
and is typically ^ 1 Gyr The full 6D information contained 
in each sample is transformed into the set of observables of- 
ten used to represent t he stream: position on the sky (A©, Bq) 
dMaiewski et al.l2003 ^. heliocentric distance d, line-of-sight ve- 
locity in the Galactic standard of rest Wgsi , and proper motions 

(/i6,/iiCOs6). 

In total 5 X 10** initial conditions are integrated, producing 
5 X 10^ points in the space of observables , which are assigned to 
a grid using the Cloud in Cell algorithm (iHocknev & EastwoodI 
[T988). Fig. |2] shows the projected density for different observ- 
ables as a function of the angular coordinate A©: P(o, A©), 
witho = {wgsr, i?0, d, /Xf,,/ii COS&}. In each panel we marginal- 
ize the density over the observed quantity o at fixed A©, that is 
P{o\Kq) = J doP{o, A0). The solid black line shows the me- 
dian of P(o|A0), and with gray bands we represent the 1 and 
2a equivalent scatter around the median. 

For comparison we have included the mean orbit of the 
LMIO model (orange dashed line), and their A^-b ody run (green 
dots). We have also added the measurem ents of Majewski et alj 
(l20()4j) (cyan stars), Correnti et al.l (12010.) (magenta triangles) 
and ICarlin et al.l (1201 2|) (red diamonds). As mentioned before, 
there are some deviations between the mean orbit and the lo- 
cation of the tidal stream as probed by the A^-body run, for 
example in the distances d of the trailing arm. 

Fig.|2]shows that the radial velocities Ugsr, distances and the 
positions in the sky Bq are well fit in our new potential, and 
as well as in the LMIO model. In test runs we found that the 
dependence of the fits on the parameter 7^ is not strong, when- 
ever this is kept within reasonable values. Of course, a value 



of Ta > 



will lead to potential that is purely oblate in the 



region probed by the stream, and therefore will not be able to 
fit the velocities of the leading arm. 

The dependence on the flattening qz, is shown in Fig. [3] for 
the leading arm (the trailing arm is rather insensitive in the re- 
gion where observations are available). We explore four differ- 
ent values of q^ = {0.7, 0.8, 0.9, 1.1} keeping Tq — 30 kpc. 
In the regions probed by the data, Aq > 200°, the effect of 
changing qz is strong on the velocities, which clearly rule out 
qz < 0.9. On the other hand, the positions on the sky disfavor 
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qz > 1. In general, we find that 0.90 < Qz 
fits to the observables in the leading arm. 



< 0.95 yield good position of the LMC, tlmc^ is 



3. OUTER HALO; THE EFFECT OF LMC 

It is very intriguing that the direction of the major axis of the 
LMIO potential approximately lies in the direction towards the 
LMC. This suggests that the LMIO potential may perhaps be 
seen as an effective field: the result of the combined potentials 
of the LMC and of a truly triaxial MW halo. 

Let us consider the various torques exerted on the (instanta- 
neous) plane of motion of Sgr First, note that since <j> w 90°, 
the principal axes of the potential of the halo are nearly aligned 
with the Galactocentric coordinate system. Consequently, we 
can simplify Eq. ^ to 
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The torque induced by the LMIO potential is simply r = r x 
9$haio/5r. Of the three components of this field, the x and z 
components are controlled by gradient of the force along the y 
direction, i.e. that of the major axis of the LMIO halo. Consider 
for instance the z component, 
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Note that in our reference system the present day position of 
the LMC is nearly on the plane a; = (See Fig. [T]). The force 
generated at r = xi + yj + zfc by a point mass A/lmc at the 
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Fig. 2. — Radial velocity Dgar, position in the sky Bq, and lieliocentric 
distance d as function of the angular distance along the stream Aq for the 
leading arm (left) and trailing arm (right) for the potential described in Fig.[T] 
The solid black line is the median orbit and the shaded regions represent 1 a 
and 2t7 equivalent dispersion. The green points are from the A'^-body simula- 
tion by LMIO, while their center of mass orbit is the orange dashed curve. 
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that generates a torque t^^'~^ — rx .Flmc^ whose z component 

is 
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Using Eqs. (O and ( fTOl i we can quantify the relative amplitude 
of the torques exerted by the triaxial halo and by the LMC on 
the orbit of Sgr at its present location 

tLmc Mlmc r^ WLMc 1 
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The mass of the LMIO halo enclosed at the present distan ce of 
Sgr is Mhaio ^ 10" Mq « AfLMC JBeslaet alj|2010h . At 
the present day r/rggr/LMC ~ 0.5, while yLMc/v ~ 10, and 
taking qi = 1.38, this implies that the expression above is of 
order unity. Additionally, since qi > I the torque generated 
by the LMC points in the same direction of that induced by the 
triaxial halo (j/lmc < 0, c.f. Fig.[T]). This means that presently 
the torque on Sgr generated by the LMC is as important as the 
one generated by the triaxial halo. 

To confirm this order of magnitude argument we perform 
new orbital integrations in a slightly modified halo model, 
which is still given by Eq. ^ but now with axis ratios qi = 1.1, 
q2 = 1.0 and q^ — 1.25, and where we have kept the orienta- 
tion (f) = 97°. We also include the potential of the LMC. To this 
end, we evolve backward and forward the orbit of the LMC, 
from its present day position (a, S ) = (5*^, 27.'^6, -69°, 52.2') 
(Piatek, Prvor& Olszewski 2008), heliocentric distance d = 
50.1 kpc (Freed man et ak ^ 001). proper motions (/i; cos b , Hb) 
= (1.96, 0.44) mas vr~^ (iPiatek Pryor & OlszewskilboOSl) and 
line-of-sight velocity Vr 



270 km s~^ (Ivan der Marel et ai] 
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Fig. 3. — Leading arm line-of-sight velocities (top), position on the sky 
(middle) and heliocentric distances (bottom) for different qz and ra = 30 
kpc. The potential for r 2> r'a is the same triaxial model as in Fig. |2] The 
dashed orange line is the mean orbit of LMIO. 
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l2002h . We then place a Hemquist sphere of mass Mlmc,o ~ 
8 X 10^° Mq and scale radius tlmco = 2 kpc along this orbit. 

Fig.|4]shows a model of the Sgr stream orbital path in which 
the LMC and our slightly revised halo are included. As before 
we show the LMIO model with the dashed orange line. It is 
interesting to note that after including the LMC in our model, 
the orbital tilt of the leading arm is well fit despite the change 
in shape of the MW halo, which is now elongated in the z- 
direction at large radii. The torque of the orbital plane is, of 
course, also felt by the trailing arm, resulting on a slight change 
in the direction of gradient of Bq for A < 100°. 

This analysis shows that the perturbations of the LMC on 
the orbit of the Sgr stream are non-negligible. In a more re- 
aUstic scenario including dynamical friction, the LMC might 
have been even more massive than at present day, and as a con- 
sequence, its role in shaping the orbit of Sgr even more im- 
portant. However, some caution is necessary before drawing 
strong conclusions about the dynamics of the stream 3-4 Gyr 
ago. For example, if the LMC is in its first infall, in which case 
the closest encounter with the Sgr dwarf galaxy is currently tak- 
ing place. However, during the last ^ 2 Gyr its presence could 
have significantly affected older wraps of the Sgr stream. 

We explore in Fig. |5]how the differences between the var- 
ious models may become apparent for older portions of the 
stream. Here we show the first (solid line) and second (dot- 
ted line) wraps of the leading (left) and trailing (right) arms, for 
the three different models discussed so far: black is our fiducial 
model from Section|2] green is the model that includes the LMC 
and red is the LMIO triaxial model. We have included observa- 
tions of different stellar tracers: RR Lyraes dlve zic et al. 2000', 
IVivas. Zinn & Gallarll2005l:lPrior. Da Costa & Kelled2009l) . Car- 
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bon Giants ftbataet al.ll200ll), RGB stars dPohm-Palmer et alj 
2001; Starkenburg et al. 2009: ICorrenti et al.ll20I0l) , M Giants 
(Mai ewski et al.l2004 ) and Red Horizontal Branch stars dShi et alJ 
2012). It should be noted that (Shi etal. 2012) preselect their 
sample according to the LMIO model. 

We show also the positions in the sky for the bright (orange 
filled squares) and the faint (orange op en squares) streams i n 
the Southern Galactic hemisphere from iKoposov et aTl (|20I2|) . 
Whereas the association to the trailing arm is clear for the 
brighter portion of the stream, the faint parallel stream could 
perhaps be an older wrap from either trailing or leading arm. 
More information, especially kinematic is necessary to disen- 
tangle the various contributions of Sgr in this region of the sky, 
and these might also help constrain further the shape of the dark 
halo of the Milky Way. It should be born in mind that although 
the differences between older wraps amongst the various mod- 
els are larger than for younger streams, the predictions for their 
properties are clearly much more uncertain. 

4. CONCLUSIONS 

In this paper we have presented a new model for the MW 
dark matter halo that fits the observations of the Sgr stream, 
both the radial velocities as well as the orbital tilt of the lead- 
ing arm. The dark halo potential is axisymmetric and flattened 
towards the disk plane for r < 10 kpc, with q^ = -9, and 
asymptotically approaches the lLaw & Maiewskil(l2010h triaxial 
model at larger radii. A gratifying property of this potential is 
that its inner oblate shape and orientation account for the pres- 
ence of the Galactic disk and ensure its stability. 

The triaxial part of this potential, however is not entirely 
consistent with expectations from the ACDM model. Its odd 
(nearly oblate) configuration can be changed, and brought to a 
more cosmologically plausible shape, if the gravitational field 
generated by the Large Magellanic Cloud is taken into account. 
The integration of orbits in a composite potential including the 
LMC and an outer triaxial halo with qi = 1.10, (72 = 1-00 and 
qs — 1.25 (that is as before, oblate in the inner regions), is also 
found to reproduce well the properties of the Sgr streams in the 
region where these have been constrained observationally. 

The conclusions drawn in this work are based on heuristic 
searches of the high-dimensional parameter space that charac- 
terizes the gravitational potential of the Milky Way and that 
of the LMC. By no means do they represent best fit models 
in a statistical sense. Therefore, the predictions made cannot 
be considered exclusive or definitive, but serve to guide where 
future observations could focus to distinguish between various 
models. Notwithstanding these caveats, we have been able to 
demonstrate that the dynamics of the Sgr streams can be under- 
stood in a context that is consistent with expectations from the 
ACDM model. 
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Fig. 4. — Radial velocity i)gsr, position in the slcy Bq, heliocentric distance 
d and proper motions fii cos b, /if, as a function of the angular distance along 
the stream Aq for the leading arm (left) and trailing arm (right). The poten- 
tial used includes the LMC as well as that for the halo, which has the form 
described in Eqs. J2)-(6), i.e. it is oblate in the center with q^ = 0.9 and 
ra = 30 kpc, but with axis ratios qi = 1.1, 52 = 1-0 and 53 = 1.25. 
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